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Abstract 

The gyraton model describes a gravitational field of an object moving with the velocity of light 
which has finite energy and spin distributed during some finite time interval L. A gyraton may be 
considered as a classical toy model for a quantum wave packet of high-energy particles with spin. 
In this paper we study a head-on collision of two gyratons and black hole formation in this process. 
The goal of this study is to understand the role of the gravitational spin-spin interaction in the 
process of mini black hole formation in particle collisions. To simplify the problem we consider 
several gyraton models with special profiles of the energy and spin density distribution. For these 
models we study the apparent horizon (AH) formation on the future edge of a spacetime region 
before interaction. We demonstrate that the AH forms only if the energy duration and the spin are 
smaller than some critical values, while the length of the spin distribution should be at least of the 
order of the system gravitational radius. We also study gravitational spin-spin interaction in the 
head-on collision of two gyratons under the assumption that the values of gyraton spins are small. 
We demonstrate that the metric in the interaction region for such gyratons depends on the relative 
helicities of incoming gyratons, and the collision of gyratons with oppositely directed spins allows 
the AH formation in a larger parameter region than in the collision of the gyratons with the same 
direction of spins. Some applications of the obtained results to the mini black hole production at 
the Large Hadron Collider in TeV gravity scenarios are briefly discussed. 

PACS numbers: 04.70.Bw, 04.30.Nk, 04.25.Nx, 04.50.+h 



I. INTRODUCTION 



The black hole formation in high-energy particle collisions is an important issue especially 
in the context of TeV gravity scenarios l|, Ij, |3|. In the theory with large extra dimensions, 
the Planck energy could be of the order of TeV and collisions of particles with the center- 
of-mass energy greater than the Planckian one will occur in future accelerators such as the 
Large Hadron Collider (LHC) at CERN [4] . A detailed study of mini-black-hole production, 
especially at the threshold of this effect, would require the complete theory of quantum 
gravity. However, if the mass of a created black hole is much larger than the Planck mass, 
one can use a semiclassical approximation to describe such processes. In this approximation, 
the black hole formation in particle collision and its subsequent decay in the process of the 
Hawking evaporation are studied in the framework of the (semi) classical general relativity. 
The apparent horizon (AH) is a useful tool for estimation of black hole production rate in 
this approximation, because the existence of an AH is a sufficient condition for the black 
hole formation. 

The first work along this line was done by Eardley and Giddings [sl in the four- dimensional 
case. Since the high-energy particles are relativistic, they used the Aichelburg-Sexl (AS) 
metric to describe the gravitational field of such particles before their collision. The 
AS metric can be obtained by boosting a Schwarzschild black hole to the speed of light and 
keeping the energy p of the boosted black hole fixed. The gravitational field of the AS particle 
is a shock and it is localized on the null plane {u = for one of the particles and v = for 
the other one). One of the null generators on each of the null planes represents a particle 
trajectory, while its gravitational field is distributed in the transverse plane orthogonal to 
the direction of motion. Two AS particles do not interact before the collision and the metric 
outside of the interacting region is known explicitly. Eardley and Giddings analytically 
studied the AH on some slice [u = > v and v = > u) and derived the maximal impact 
parameter femax for the AH formation. The quantity o"ah = "^^^max gives the lower bound on 
the cross section of the black hole production. 

The results of js] were generalized by one of us and Nambu [t^ for the mini-black-hole 
formation in the higher-dimensional spacetimes. In this work the numerical calculations were 
used. Later, one of us and Rychkov [s] improved these results by studying the AH on the 
futuremost slice that can be adopted without the knowledge of the interacting region (i.e.. 
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u = < V and v = < u). Using this approach the stronger lower bound (JAH/7r[rft,(2p)]^ ~ 
2-3 on the cross section of the black hole production in the collision of AS particles was 
obtained for spacetime dimensionality D = 5-11. 

Certainly the model of colliding AS particles is oversimplified. The AS particles are 
assumed to be neutral and spinless. In reality all the known elementary particles have spin 
and most of them have either electric or color charge as well. 

Charged particles have additional charge-charge interaction besides the gravitational in- 
teraction. Moreover their gravitational interaction can also be modified. The latter effect 
was discussed to some extent by one of us and Mann ^ . In that paper a boosted Reissner- 
Nordstrom metric was used as the model of an ultrarelativistic charged particle and the 
head-on collision in a spacetime with an arbitrary number of dimensions was studied. The 
results obtained in that paper indicate that the charge makes the AH formation more diffi- 
cult. It was also argued that the effects of the quantum electrodynamics could change the 
results. The results of {q] were later used by Gingrich 10 1 who reconsidered the black hole 
production rate at the LHC. 

In the quantum mechanical description, the colliding particles are characterized by wave 



packets which have finite duration in time 



ll| . To take into account this effect as well as 



to include spin-spin interaction, in this paper we study head-on collisions of two gyratons. 

The gyraton model was proposed in [l2|. The motivation of this paper was to find 
the gravitational field generated by a beam pulse of spinning radiation with a finite time 
duration, which is propagating at the speed of light. In the gyraton model, the metric 
outside of the source satisfies the vacuum Einstein equations and the gravitational field is 
distributed in the plane transverse to the direction of motion. Unlike an AS particle, the 
gravitational field of a gyraton is not a shock wave but has the finite duration in time. 
A gyraton may have spin which manifests itself in the dragging-into-rotation effect. The 
AS particle metrics can be obtained from the gyraton solutions if the duration is taken 
infinitely small and the spin vanishes. General properties of gyraton metrics were studied in 
detail i n |l3l] . Electrically charged gyratons and gyratons in the supergravity were discussed 
in 13], respectively. Gyraton solutions can be also generalized to the case when the 
spacetime is asymptotically anti-de Sitter [l3| . 

Colliding gyratons which we consider in this paper differ from the AS particles both 
by the presence of spin and the finite duration in time. Let us discuss briefiy what kind 
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of new features one can expect. First natural question is: Can one include spin effects 
in the interaction between highly nonrelativistic particles by boosting the Kerr metric? 



dimensional spacetime and in 



Such boosted Kerr black hole solutions were considered, e.g. in [18|, ll9|, [20I, in the four 



2l| in higher dimensions. The main problem in this approach 



is the following. In order to have a well-defined limit for the boosted metric, one needs to 
keep the energy p of the system fixed, so that the mass of the black hole M = 7~^p must 
vanish when the 7-factor infinitely grows. If one assumes that the spin s remains finite in this 
limit, the rotation parameter a = s/M infinitely grows, so that the metric describes a naked 
singularity. The radius of the ring singularity is of the order of a and also infinitely grows. 
The latter problem can be avoided by assuming that the rotation parameter a remains finite 
in the infinite boost limit, as it was done in the above references. Although finite results 
different from the AS particle can be obtained by this procedure, fixing a means that the 
spin s = aM of the boosted object vanishes. Furthermore, in this limit we have an object of 
typical size a, which does not satisfy the requirement that we would like to have a point-like 
object. Thus a boosted Kerr black hole does not provide one with a suitable model for an 
ultrarelativistic particle with spin, e.g. for a photon. In the gyraton model the spin is easily 
included. 

There is another aspect of the high-energy particle collision which the gyraton model may 
help to understand better. Recently the validity of an AS particle as the model of a high- 
energy elementary particle was questioned in |22:]. In this model the curvature invariants 
at the moment of the collision of the two planes, representing the colliding particles, are 
infinite, so that formally higher-curvature quantum corrections may be important. This 
problem was discussed in llj . It was argued that the quantum effects, such as the finite size 
and finite duration in time of the incoming wave packets, can help to solve this problem. 
The quantum-to-classical transition in the description of the mini-black-hole formation in 
the particle collision is an interesting open question. We do not address it in our paper, 
but instead we use the gyraton model in order to discuss how the finiteness of the duration 
in time of the colliding objects modifies the results of the AS approximation. In such an 
approach, the gyratons might be considered as an effective model for the quantum wave 
packets. 

With these motivations we study the AH formation in the head-on collision of gyratons. 
For simplicity we consider the four- dimensional case. The paper is organized as follows. In 
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the next section, we introduce several gyraton models: a gyraton without spin, its AS limit, 
and spinning gyratons. As for spinning gyratons, we consider two types depending on relative 
locations of the energy and spin profiles. Then we set up five cases of head-on collisions of 
gyratons. In Sec. Ill, we derive an AH equation on the future edge of the spacetime region 
before interaction, u = < v and v = < u. We present the numerical results of this study 
in Sec. IV. The conditions on the spin and on the energy and spin durations for the AH 
formation are obtained for each of the collision cases. Then in Sec. V we focus our attention 
on the study of the spin-spin interaction. In general, this is a complicated problem, since it 
requires the knowledge of the metric in the region of interaction. To obtain it, one needs to 
solve nonlinear Einstein equations. We simplify the problem by assuming that the spins of 
the interacting objects are small and solve the equations by using a method of perturbation. 
Then we study again the AH formation on the new slice that is the future edge of the solved 
region. In the adopted approximation we obtain spin-spin interaction corrections to the 
mini-black-hole production. Sec. VI contains a summary of the results and a discussion of 
their possible applications for the study of mini-black-hole production at the LHC 

II. SYSTEM SETUP 

In this section, we set up the problem of the head-on collisions of two gyratons. We first 
review the gyraton model in a four-dimensional spacetime in Sec. II A. The gyraton has 
the total energy p and the spin J which are distributed in time. Their time profiles are 
characterized by two functions. We introduce four kinds of gyratons by specifying these two 
functions. Next in Sec. II B, we introduce a system of null geodesic coordinates, which is 
necessary for specifying the slice to study the AH existence. It is also useful for clarifying 
the gravitational field of gyratons. In Sec. II C, we set up five cases of head-on collisions of 
two gyratons, using the introduced four gyratons. 

A. The gyraton model 

The gyraton model was proposed in [l^. In that paper, the gravitational field generated 
by a beam pulse of spinning radiation was first studied in the weak field approximation 
and then the exact solutions of the Einstein equations were obtained, which reduce to the 
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approximate solution at far distance from the source. These solutions were obtained in any 
number of spacetime dimensions. In particular, a four- dimensional gyraton has the metric 

ds^ = —dudv + df^ + f^(i0^ — AG (2p\ogfXp{u)du — Jxj{u)d(j)) du. (1) 

This metric represents a spacetime in which a segment-shaped source located at f = with 
the total energy p and the spin J is propagating at the speed of light along u = const. 
The existence of the term d(f)du in Eq. ([T]) indicates presence of a dragging-into-rotation 
effect generated by the spin source. The functions Xp{^) ^"^^ Xj{^) ^i-re arbitrary functions 
satisfying the normalization conditions 

j Xp{u)du = j Xj{u)du = 1. (2) 

They represent the energy and spin density profiles, respectively. 

Introducing the new coordinates x := (f cos0, f sin0), non-zero components of the 
energy-momentum tensor of the gyraton are calculated as 

Tnu = VXp{u)5{x) + T,GJ\]{u)5\x), (3) 

Tua = ^Xjiu)eabdb5{x). (4) 

These formulas show that the source has an infinitely narrow shape. For a realistic beam 
pulse, the source should have a finite radius f = f s and the metric will take a different form 
from Eq. ([1]) for f < fg. Therefore Eq. ([1]) is considered to be the metric outside of the 
source f > fg. In this paper, we do not take account of this finiteness (in space) of the beam 
pulse for simplicity and adopt Eq. ([T]) for arbitrary values of f (i.e., fg = 0). 

Hereafter we adopt the gravitational radius of 2p, i.e. rh{2p) = AGp, as the unit of the 
length. In this length unit, the gyraton metric is represented as 

ds'^ = —dudv + df^ + f'^dcj)^ — 2 \ogfXp{u)dv? + 2jxj{u)d(j)du. (5) 

Here, j is a dimensionless quantity defined by j := J/2prh{2p) and we use j as a parameter 
to specify the spin of the gyraton. 

The gyraton model is specified by determining the functions Xp(^) ^iiid Xj{^)- The inter- 
action between two gyratons with arbitrary profiles Xp{^) ^"^^ Xj{^) is a quite complicated 
problem which, in the general case, requires numerical calculations. Hence it is natural to 
consider first the simplest profiles for which the null geodesic coordinates can be studied 
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FIG. 1: The energy and spin density profiles, Xp(^) and Xji^), for the p-, AS-, a-, and b- gyratons, 
respectively. The gray arrows show the directions of propagation. 

analytically. For this reason, in this paper we consider four types of gyratons whose energy 
and spin profiles are as shown in Fig. [H We will explain them one by one in the following. 
For convenience, we introduce the following function 

i}iu,L) = jieiu)-eiu-L)) , (6) 

where d{u) is the Heaviside step function. Its integral over u is 1, and in the limit L — > it 
gives a 5-function. 



1. p-gyraton 

The first one is a gyraton without spin with energy duration L. For this model, we adopt 

Xp{u) =^{u,L), 
Xj{u) = 0. 

This model is useful for studying the effect of the energy duration on the AH formation. 
We simply call it a "spinless gyraton" or a "p-gyraton" because it has only one parameter, 
energy p. 



2. AS-gyraton 

The second one is an Aichelburg-Sexl (AS) particle [6j with 

Xp{u) = 6{u), 



(8) 

Xj{u) = 0. 



This is the hmit L — > of the p-gyraton. Hereafter we call it an "AS-gyraton" for short. 

3. a-gyraton and b-gyraton 

The remaining two gyratons, which are referred to as an "a-gyraton" and a "6-gyraton", 
have nonzero spin. We adopt the following functions of Xp('u) and Xji''^) fo^ ^-^d ^ 
gyratons: 



In these two models, L represents the spin duration and the energy has zero duration. We 
call them an a-gyraton and a 6-gyraton, respectively, because for the a-gyraton, the spin 
source comes after the energy source, while for the 6-gyraton, the spin source comes before 
the energy source. These two models are useful for studying the effect of the spin and its 
density duration on the AH formation. By comparing the results of a- and b- gyratons, 
we can understand to what extent the obtained results depend on relative positions of the 
energy and spin density profiles. Each of the gyratons is reduced to an AS-gyraton if we 
take j = 0. 

Readers might wonder why we do not adopt Xp('u) = Xji"^) — ''^{'^^ L) for spinning gyra- 
tons. This is because of a technical problem. In the next subsection, we derive a coordinate 
system based on the null geodesic congruences. This coordinate system can be analytically 
derived for a- and h- gyratons, but not for Xp(^) = Xj{^) = '^{u, L). 




(10) 
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B. Null geodesic coordinates 



In this subsection, we introduce null geodesic coordinates, which are very useful for 
specifying the slice on which we study the AH. We introduce new coordinates {u, v, r, 0) by 

u = u, 

V = v + F{u,r), ^^^^ 

f = G{u,r), 

= + H{u, r). 

We assume that the two coordinate systems coincide for u = u < and hence F = H = 
and G = r for u < 0. For m > 0, we require a line v,r,(f) = const, to be a null geodesic and 
the coordinate u to be its affine parameter. This requirement is realized if and only if the 
following relations are satisfied: 

Hu = (12) 

F„ = Gl - 2xAu) logG - ^-gf^, (13) 

= 2G,„G,,,. (14) 

These relations determine F, G, and H. In terms of these functions the metric takes the 
form 

ds^ = -dudv + G^dr'^ + G2(d0 + H^rdrf. (15) 
Eliminating F from Eqs. (fT3l) and (JT4l) . we find 

^ _ Xp{u) 3\]{u) 

<^,uu ^ 1 ^ • 



Once this equation is solved, one can find H by solving Eq. (1121) and determine all the 
coefficients in the metric (ITSl) . 

1. p-gyraton and AS-gyraton 

For a spinless gyraton, using Eq. ([7]), we find 

, |G(n,r), (0<«<L), 
G'(M,r)=<^ (17) 

G„(L,r)(^z-L) + G'(L,r), (L < m). 
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FIG. 2: The left plot shows the coordinate singularity in the (u, r) coordinates where G = for 
a p-gyraton with L = 1. The lines r = const, are null geodesies and they hit the coordinate 
singularity. The right plot shows the light rays in the (n, f) coordinates. The light rays with an 
identical r value bend due to an attractive force and focus to one point on a symmetry axis, which 
corresponds to G = 0. 

where 

G(m, r) := r exp [- [erf" ^ {y)] ^) , (18) 



Here, the function erf~^(7/) denotes the inverse function of the error function: x = eri~^(y) 
is equivalent to 

2 

y = eTi{x):= —exp{-f)dt. (20) 



In the limit L —>■ 0, the function G reduces to 

G{u,r) = r-u/r, (0 < m), (21) 



and the metric f|T5|) coincides with the AS-gyraton in the null geodesic coordinates |8l. |22[| . 
We should point out that there is a coordinate singularity at 

\^^r, iO<u<L), 
u = < ^ _ (22) 

[L-G(L,r)/G,(L,r), (L < m), 

where G = 0. The shape of the singularity is shown by a solid line at the left plot of Fig. [21 
In order to understand the meaning of this singularity, it is useful to consider null geodesies 
v,r,(j) = const. These null geodesies plunge into the coordinate singularity. Let us go back 
to the original {u, v, f, 0) coordinates. The trajectories of the light rays in the coordinates 
{u, f) are shown in the right plot of Fig. [2l Because of the gravitational effect of the gyraton 

10 



u 



u 



FIG. 3: The left plot shows the coordinate singularities in the (u, r) coordinates where G = 
(dark gray lines) and G^r = (light gray line) for an a-gyraton with L = 1. The lines r = const, 
are null geodesies and they hit one of the two coordinate singularities. The right plot shows the 
propagation of light rays in the (n, f) coordinates. Light rays with a large r value focus to one 
point, which corresponds to G = 0. On the other hand, light rays with a small r value bend 
outward due to a repulsive force around the center. Then two neighboring light rays cross each 
other, which corresponds to G^r = 0. 

energy, the proper circumference of a congruence of light rays with an identical r value 
becomes small as v is increased and eventually becomes zero. This is where the congruence 
hits the coordinate singularity in the (m, r, cf)) coordinates. Thus, the coordinate singularity 
corresponds to the symmetry axis and therefore we call it the focusing singularity. 

2. a-gyraton 

Now we turn to the spinning a-gyraton. Using Eq. ([9]), Eqs. (|T6|) and ([T2|) are solved as 




(0 < M < L), 



(23) 




(0 < M < L), 



(24) 



where 




(25) 




(26) 
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FIG. 4: The same as Fig. [3] but for a 6-gyraton. 



In this case, there are two coordinate singularities. One is the singularity at 

1 + jVL2 - r'^/L 



.^ /^'!-r^g,. . (IS.). (27, 



where G = 0, and the other is at 



u = (28) 

L-G,(L,r)/G,,(L,r), (L < m), 

where G = 0. The two coordinate singularities are shown in the left plot of Fig. [3l A 
light ray v,r,(j) = const, plunges into one of the two singularities. The propagation of light 
rays in the {u, f) coordinates is shown in the right plot of Fig. [3l Because there is the 
energy distribution at m = 0, the light rays bend quickly there. Then the light rays with 
a sufficiently large r value focus to one point and this is the focusing singularity G = 0. 
On the other hand, around the center, the gravitational field generated by the spin source 
is repulsive and the light rays of a small r value bend outward. Because of this effect, 
the neighboring light rays cross with each other and this is where the congruence hits the 
coordinate singularity G^r = 0. Therefore we call it the crossing singularity. 



3. b-gyraton 

Finally we obtain the formulas for G and H of a spinning 6-gyraton. They are found by 
solving Eqs. ( |T2ll and (|T6ll using Eq. ( ITOl) . The result is 

G{u,r), (0<M<L), 
12 



H(u,r), (0<u<L), 
H{u,r) = { ' - - (30) 

H{L,r), {L<u), 



where 



^K-)--yi+P7^; (31) 

H(u,r) := — arctan . (32) 



Similarly to the a-gyraton, there are the focusing singularity at 



r2(rVL+l) , , 
^= ryi-f ' (^^^)' (33) 

and the crossing singularity at 

AL/\j\)r^, (0<w<L), 
M = { (34) 

L-G,(L,r)/G,,(L,r), (L < m). 

The shape of the two singularities in {u, r) coordinates and the propagation of light rays in 
[u, f) coordinates are shown in Fig. HI 



C. Gyraton collisions 

Consider two gyratons and assume that each of them belongs to one of the four types 
described above. We obtain several different configurations for the collisions of these gyra- 
tons. To illustrate main features of these collisions, in this subsection we set up five cases of 
the head-on collision of two gyratons which are of the most interest. The incoming gyratons 
are referred to as gyraton 1 and gyraton 2. Let us divide a spacetime for the two-gyraton 
system into four regions: 

Region I: (u < 0,v < 0), 

Region II: (u>0,v<0), 

^ ~ ^ ^ (35) 
Region III: {u<0,v>0), 

Region IV: {u>0,v>0). 

Because the gyratons propagate at the speed of light, they do not interact with each other 

before the collision. Thus we can use the metric of the gyraton 1 in the regions I and II and 
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FIG. 5: The five cases of gyraton collision that we study in this paper. In case (0), two identical 
p-gyratons collide. In cases (la) and (lb), a- and b- gyratons with spin j collide with an AS 
particle, respectively. In cases (2a) and {2b), two a- and b- gyratons with spins j and aj (cr = ±1) 
collide, respectively. 

the metric of the gyraton 2 in the regions I and III (by changing u and v). Then, the metric 
of the system is given as 



—dudv + dr"^ + r'^dcj)^, 

-dudv + G^P{u,rfdr^ + G^^\u,rf{d<i) + H^r\u,r)drf 
-dudv + rfdr"^ + G^^K^, r)'^{d<P + Hf{v, r)dr f 



(Region I), 
(Region II), (36) 
(Region III). 



The metric of the region IV is unknown, because the interaction between the two gyratons 
determines its structure through the Einstein equations.^ 

In the previous subsection, we introduced four gyraton models, i.e., a p-gyraton, an AS- 
gyraton, an a-gyraton and a fe-gyraton. Using these models, we will consider five cases of 
collision. The first one, which we call the case (0), is the collision of two identical spinless 
p-gyratons. For both G*^^-* and G^'^\ we use the formula of G for the p-gyraton f[T7l) - ffT^ . 
The energy p determins the scale, so that the only essential parameter which specifies the 
system is the energy duration L. 

In the next two cases (la) and (Ife), we consider collisions of a spinning gyraton (a- and 
h- gyraton, respectively) with an AS-gyraton. In both cases, we assume that incoming 
gyratons have the same energy, and only a gyraton 1 has the spin j. These are interpreted 
as collisions of a spinning particle and a particle without spin. For G^^-* and H'^^\ we use 
the formulas fl25]) - fl2^ of G and H for a-gyraton in the (la) case, and use the formulas 
([29])-([32]) of G and H for ?^gyraton in the (16) case. For G^'^\ the formula ^ oi G for 



^ In Sec. IV, we will solve the Einstein equation in a part of the region IV in a specific case, when spins 
are small. 
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AS-gyraton is used in both cases. The essential parameters which specify the system are 
the spin j and its duration L for a gyraton 1. 

In the remaining two cases (2 a) and (26), we study colhsions of two a-gyratons and two 
b- gyratons, respectively. These are interpreted as collisions of two spinning particles. In 
these cases, the incoming gyratons are assumed to have the same energy p and the same 
spin duration L. As for the spin values, we assume that the gyraton 1 has the spin j and the 
gyraton 2 has the spin aj, where a = ±1. Therefore two spins have the same absolute value 
|j| and have either the same sign or different signs. For the choice a = +1 two spins have 
the same direction (i.e. helicities have opposite signs), and for the choice a = — 1 two spins 
have opposite directions (i.e. helicities have the same sign). For G'(i) and H^^\ we use the 
formulas fl2^ - fl2Bl) of G and H for a-gyraton in the (2a) case, and use the formulas fl2^ - fl52]) 
of G and H for 6-gyraton in the (26) case. For G*^^-* and H^'^\ we use the formulas of G 
and H for a- and b- gyratons with j replaced by aj in the (2 a) and (26) cases, respectively. 
The essential parameters which specify the system are the spin j of gyraton 1, the relative 
directions of two spins a, and the spin duration L of each incoming gyraton. In the study 
of Sees. III-IV, the condition for AH formation in the slice u = < v and v = < u will 
turn out to be independent of a, and hence the essential parameters are reduced to j and 
L. The sign of a will become important in the study of the spin-spin interaction in Sec. V. 

All the five cases are schematically illustrated in Fig. [51 

III. FINDING THE APPARENT HORIZON 

The apparent horizon (AH) S is a compact two-dimensional spacelike surface for which 
the family of outgoing null rays emitted orthogonally to S has zero expansion. We study 
the AH on the slice u > = v and v > = u. Figure shows a schematic picture of the 
AH for colliding gyratons. The AH consists of two parts: 



These parts are connected at u = v = 0. Each surface has the other end at the focusing 
singularity. Because the focusing singularity is just one point for the same r value, the 
surface given by Eq. ( 1371) is a two-dimensional closed spacehke surface. 





(t) > = u). 



(37) 
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FIG. 6: The schematic picture of the AH in the shce u >0 = v and v > = u. 



Because the AH equations and the outer boundary conditions for h^-^^ and h^"^^ have the 
same form, we only consider h{r) := h^^\r) and denote G = G^^^ and H = H^^\ The metric 
in the neighborhood of i; = < m is given by 



ds^ = -dudv + G^^dr^ + G\d(j) + H^rdrY 



(38) 



{h{R),R,$). The local lightcone with the apex at this 

(39) 



Let us consider a point (m, r, (j)] 
point is 

{u - h{R))v = Glir - Rf + G\<P - <P) + H^r{r - i?)]'. 

We find the envelope of the local lightcones by taking the derivative of Eq. (I39ll with respect 
to R and The tangent vector of the null geodesic congruence in the {u, v, r, 0) coordinate 
is 

k^=(rS^.2.h.-^hr]. (40) 



_^ 2 -^h 

Now we calculate the expansion. The induced metric on f = const surface is given by 
d'y'^ = G^dr"^ + G'^{d(f) + H^rdr)"^, and its determinant 7 is = GrG. Let us consider a 
rectangular coordinate domain with apices at the four points Pi, P2, P3 and P4: 



(41) 





(P+ 


,^+) 




(P+ 


,^-) 


< 






(R- 


,^+) 


.^4: 


(R- 
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Here R±^ R± Ar/2 and = ± A0/2. The proper area AA{0) of this domain is 



The null geodesies passing through the apices are 

Pi 



(42) 



(r, 0) 



(43) 



(i?+ + k'{R+)X, ^+ + k't'{R+)X), 
{R+ + k'{R+)X, ^_ + k'l'{R+)X), 
{R_ + k'-{R_)X, + k'^{R_)X), 
+ r(7?_)A,^>_ + k'^{R-)X), 

where A is an affine parameter. In what follows, we keep terms up to first order in A. 
The coordinate shape of the domain surrounded by the four apices after evolution is a 
parallelogram as indicated by the vectors 



P' 

P' 

P' 
-'4 



P^P[^PiP^ = (O,A0), 



P^P[ = PiP^ = ((1 + drk''X)Ar, drk^XAr). 



(44) 
(45) 



The coordinate area of the domain is (1 + drk'^X)ArA(f) and the proper area of the domain 
is 



AA(A) = {G,rG)\ 
= A^(0) 



{h{R)+k^{R)X,R+k^ {R)X) 



;i + <9,rA) ArA0 



T , / G ruk^ + G rrk^ , G uk^ + G rk^ , , r\ \ 

+ i — g; — + — G — ^^^^)\ 



Hence, the condition dAA/dX — implies 



{3/2)G^ruh,r + G^rr / '^)G ,r + G ■ 



G 



+ 



G 



0. 



(46) 



(47) 



This equation determines the AH surface S. 

Let us discuss now the boundary conditions at the outer boundary r = Tmax- By the 
continuity of the surface, the AH should cross the coordinate singularity at r = r^ax' 



G'(/l(rniax),rmax) = 0. 



(48) 



The continuity of k^ also must be imposed, because the surface has a delta-functional ex- 
pansion if k'^ is discontinuous. Going back to the {u, v, f, 4>) coordinate, the continuity can 
be imposed as k^' = by the axisymmetry. This is equivalent to 



^,r(^max) — —'^G^r/G^,, 



(49) 
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The other condition for the continuity of k'^ is that k'^ should take a finite value. But this is 
automatically satisfied since the condition = at the focusing singularity implies that 
k^ = 0. 

Now we turn to the boundary conditions at the inner boundary r = rmin- The inner 
boundary conditions depend on both h^^^ and h^'^\ By continuity of the surface, both sides 
of the AH should cross u = v = Oatr = rmin, and thus 

h^'\r^in) = h'^'\r^in) = 0. (50) 

Also the null tangent vectors k^-^^'^ and k^'^^^ of two surfaces should be parallel at r = rmin 
so that there is no delta-functional expansion, k^'^^ and k^^^" are given by 

k^''" = (2, j , = (^^, 2, h^?, j , (51) 

and k^^^ \\ k^"^^ is equivalent to 

/i.l'Hr^in)/i?(r^in)=4. (52) 

The numerical procedure for defining the AH is straightforward. First we choose some 
value of Tmax and solve h{r) with the outer boundary conditions (jUj) and fH9l) . When h{r) 
becomes zero at r = r^^^, we check whether the inner boundary condition (1521) is satisfied. 
Iterating these steps for various values of rmax, we determine whether the AH exists and find 
its location. 

Note that H does not appear in the AH equation and the boundary conditions. This 
means that the dragging-into-rotation effect causes a change in the shear but not in the 
expansion. Thus on the slice we have adopted, the condition for the AH formation does 
not depend on the sign of j in the cases (la) and (lb). In cases (2a) and {2b), it does not 
depend also on the relative directions of two spins a. Therefore in the next section j is 
assumed to be positive without loss of generality and we do not specify the value of a when 
the numerical results are shown. 

IV. NUMERICAL RESULTS 

In this section, we present the numerical results of the AH studies. The results for case 
(0), cases (la) and (16), and cases (2a) and {2b) are provided in Sees. IV A, IV B, and IV 
C, respectively. 
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FIG. 7: The top views of the AH in case (0) for L = 0.5, 1.0, and 1.423. For all < L < 1.423, we 
find only one solution. At L = 1.423, the AH almost touches the energy source at the symmetry 
axis. No AH exists for L > 1.424. 

A. Collision of gyratons without spin 

We begin with case (0), the coUision of two identical spinless p-gyratons with energy 
duration L. Figure [7] shows top views of the AH for L = 0.5, 1.0 and 1.423. For < L < 
1.423, we found only one solution, and therefore there is no inner boundary of the trapped 
region. For L = 1.423, the AH intersects the focusing singularity at m ^ L and almost 
touches the source of the energy which distributes at < m < L on the symmetry axis. 
For L > 1.424, we found no solution. Thus, on the slice we studied, the condition of AH 
formation is given by L < 1.423 in the length unit rh{2p) = 1. 

B. Collision of a spinning gyraton with an AS particle 

Next we show the results of cases (la) and (16), i.e., the collisions of spinning a- and b- 
gyratons with an AS particle. 

Figure [8] shows the AHs in case (la) for parameters L = 0.1 and j = 0.05,0.10, and 
0.1466. We found two solutions to the AH equation, which correspond to the AH and the 
inner boundary of the trapped region. As the value of j is increased for a fixed value of L, 
the trapped region grows smaller and the two solutions coincide at some value of j = jcrit- 
The trapped region vanishes for j > jcrit- The similar phenomena was observed also in case 
(16). In Fig. [9l the AH shape in case (16) is shown for the same parameters as those in 
Fig. [HI Again, there are two solutions and they degenerate at some critical value j = jcrit- 

As it has been found above, the spin makes the formation of the AH more "difficult". 
This is because the gravitational field generated by the spin source is repulsive as we pointed 
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FIG. 8: The top views of the AH in case (la) for L = 0.1 and j = 0.05,0.10, and jcrit = 0.1466. 
For all < J < 1.466, there are the AH (solid lines) and the inner boundary of the trapped region 
(dashed lines). The trapped region shrinks as j is increased and no AH exists for j > 0.1467. 
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FIG. 9: The same as Fig. [8] but for case (1&). For L = 0.1, the critical value of AH formation jcrit 
is almost the same as that in case (la). 



out in Sec. II. As the value of j is increased, the repulsive force surpasses the attractive 
force generated by the energy source and causes the extinction of the AH. 

We studied the value of jcrit as a function of L, i.e., jcrit(-^)- In Fig- [lOl the critical lines 
for the AH formation in the (L,j)-plane are shown for both cases. The AH formation is 
allowed under the critical line. For L < 1, the two critical lines agree well and go to zero in 
the limit L ^ 0. Both critical lines have the peak j ~ 0.24 at L ^ 0.6. For L > 1, the two 
critical lines show different behaviors. In case (la), jcrit{L) decreases and becomes zero at 
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FIG. 10: The critical line on (L,j)-plane for AH formation in the cases (la) and {lb). The AH 
formation is allowed under the critical line. The two critical lines almost coincide for L < 1 and go 
to zero in the limit L ^ 0. The critical line of the (la) case intersects the L-axis at L = e, while 
that of the (16) case becomes exponentially close to the L-axis as L is increased. 

L = e := exp(l). On the other hand, it decays (almost) exponentially but never becomes 
zero in case {lb). 

Let us discuss the reasoning for these results. The reason why jcrit goes to zero in the 
limit L — s> is as follows. As stated above, the extinction of the AH is caused by the 
repulsive force generated by the spin source. Thus, it is useful to introduce a radius req 
where the attractive force due to the energy and the repulsive force due to the spin balance 
one another. For this purpose, let us recall Figs. [3] and H] that show the propagation of light 
rays through the gravitational field of the gyratons. The figures indicate the existence of 
such that the rays with r > shrink and those with r < req expand in the region u > L. 
Such Teq is fouud by the equation G'^u(L,rcq) = and solved as 



eq 



L + j^/L, a-gyraton, 
^j'^/L, fc-gyraton. 



(53) 



As L is decreased, req becomes larger, which indicates that the repusive force becomes 
stronger. It is natural that req < 1 represents the condition for AH formation, and it is 
reduced to < L in the limit L — > 0. This explains the behavior of the critical lines at 
L < 1. 

At L > 1, the condition rgq < 1 does not explain the numerical results well. This 
is because the above discussion takes account only of the gravitational structure in the 
transverse direction of motion, which would be sufficient in the case L <^ 1, while the spin 
duration L plays an important role for the AH existence in the case L > 1. Let us first 
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consider case (la). Taking a limit j ^ for L > 1, the AH solution is expected to reduce 
to that for the collision of two AS particles: 



with rmin = 1 and rmax = \f&- This statement holds only for L < e. In the case L > e, 
the "solution" ( IMI) plunges into the crossing singularity. In other words, it crosses the spin 
source distributed on the symmetry axis for < u < L, on which the outer boundary 
condition cannot be imposed. Thus in case (la), the situations j = and j = 0+ are 
different. This is the reason why the critical line intersects the L 

Next we discuss case (16). In the limit j — > for L > 1, the AH solution reduces to 



because the AH never touches the spin source. Then, the condition for AH formation in the 
case j > is expected to be rmin ^ Teq, which is equivalent to < Le~^. This explains the 
exponential decay of the critical line in the (16) case. 

Note that the above interpretations, especially the ones for L > 1, strongly depend on the 
slice we have adopted. Thus there is no reason why the above discussion holds for another 
slice that is future to our slice. Hence, we should keep in mind the possibility that the 
critical line does not touch the L-axis in another slice also in the case (la). 

To summarize, for the collision of a spinning gyraton with the AS-gyraton and for the 
slice we have adopted, the condition for the AH formation is roughly expressed as L ~ 1 
and j < 0.25 in both cases. 

C. Collision of two spinning gyratons 

Finally we show the results of cases (2a) and (26), i.e., the collisions of two spinning a- 
and b- gyratons (identical up to helicities). 

Figure [TT] shows the AHs in case (2a) for parameters L = 0.1 and j = 0.05,0.10, and 
0.113. Similarly to (la) case, there are two solutions to the AH equation, which surround 



/i«(r) = h^^\r) = 2r^logr, 



(54) 




(55) 



h'^^\r) = 2r^log(r/rmm), (^mm <r< y/er^in), (56) 
where rmin = e~^/^. In contrast to the (la) case, this statement is valid for arbitrary L, 
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FIG. 11: The top views of the AH (sohd hnes) and the inner boundary of the trapped region 
(dashed hnes) in case (2a) for L = 0.1 and j = 0.05,0.10, and 0.113. No AH exists for j > 0.113. 



/-O.Ol ^ ,;-0.05 ^ j.0.0981 




FIG. 12: The top views of the AH (sohd hnes) and the inner boundary of the trapped region 
(dashed hnes) in case (26) for L = 0.1 and j = 0.01, 0.05, and 0.0981. No AH exists for j > 0.0982. 

the trapped region, and they coincide at some value of j = jcrit as the value of j is increased 
for a fixed value of L. The similar phenomena was observed also in case {2b). In Fig. [121 
the AH shape in case (26) is shown for L = 0.1 and j = 0.01,0.05, and 0.0981. Again, we 
found two solutions and their disappearance at some critical value j = jcrit- Similarly to 
cases (la) and (lb), the spin has the effect to make the AH formation more difficult. 
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FIG. 13: The critical lines on (L, j)-plane for AH formation in the cases (2a) and {2b). The two 
hnes go to zero in the hmit L — > 0. The critical lines of (2a) and {2b) cases intersect the L-axis at 
L = e and 1/e, respectively. 
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We studied the value of jcrit as a function of L. In Fig. [13], the critical line for cases (2 a) 
and {2h) in the (L, j)-plane is shown. The AH formation is allowed under each critical line. 

Let us first discuss the critical line of the (2 a) case. It goes to zero in the limit L — 
and intersects the L-axis at L = e. It has a peak j ^ 0.19 at L ~ 0.6, and this peak value 
is somewhat smaller than the peak value 0.24 of the (la) case. Hence the critical line of 
the (2a) case has the same features as that of the (la) case except that the peak value is 
smaller. For the behaviors at L ^ and L e, the same reasoning to the results of the 
(la) case holds. Compared to the (la) case, the AH formation is expected to become more 
difficult, since both gyratons have the repulsive forces around their centers in the (2a) case 
while only one gyraton has the repulsive force in the (la) case. This leads to the smaller 
peak value of jcrit (-^) in the (2 a) case. 

Now we discuss the critical line of the (26) case. It goes to zero in the limit L ^ with 
the same reason to the (lb) case. It has a peak j ~ 0.105 at L ~ 0.16, and intersects the 
L-axis at L = 1/e. The allowed region of the {2b) case is much smaller than that of the (2a) 
case. The condition of the AH formation strongly depends on the relative locations of the 
energy and spin profiles. The reason can be understood as follows. In the limit j 0, the 
AH becomes 

{2r^\og(r/y/r~) + L, {y/r~ < r < y/er^) , 
(57) 
2Tmin lo§('"/'"min) ) ('"min ^ v^'^min) ; 

where r^in is given by the equation 

^ = -^minlogrmin- (58) 

This equation has two solutions for < L < 1/e, one degenerate solution for L = 1/e, 
and no solution for L > 1/e. Thus the AH formation in the j — limit is allowed only 
for < L < 1/e. This is the reason why the allowed region is restricted toO<L<l/e 
and is much smaller than that of the case (2a). However, we should keep in mind that this 
discussion is specific to the slice we have adopted. In the case j = 0, the AH formation is 
allowed on a slice appropriately taken at the future to our slice. Hence, the allowed region 
in the (26) case is so small because of the artificial effect of the slice choice. In the next 
section, we demonstrate that this expectation is true by solving a part of the spacetime after 
the collision using the method of perturbation. 

To briefly summarize, for the collision of two spinning gyratons, the condition of the AH 
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FIG. 14: Schematic spacetime structure of the gyraton cohision in the (26) case. Region IVA can 
be solved using a perturbative method, assuming the spins are small. Then we study the AH on 
the new slice shown by a dotted line, which is the future edge of the solved region. 

formation on the slice we have adopted is roughly expressed as L ~ 1 and j < 0.2 in the 
(2a) case, while L ~ 0.15 and j < 0.1 in the (26) case. 

V. SECOND-ORDER EVOLUTION 

In a general case, finding the spacetime structure after the collision of two gyratons 
requires numerical simulations. However, in the {2b) case, we can go a little bit further 
using the method of perturbation assuming that the spins of incoming gyratons are small. 

Figure [T4l shows the schematic spacetime structure in the (26) case where the gyraton 1 
with the spin j and the gyraton 2 with the spin aj collide. (Here a = ±1, but we note that 
the solution to the Einstein equations found in the Sec. IV A can be applied to an arbitrary 
value of cr.) For a > the two spins have the same direction (i.e., the helicities of incoming 
gyratons have opposite signs), while for a < the two spins have opposite directions (i.e., 
the helicities of incoming gyratons have the same sign). The exact metrics in the regions 
I, HA, IIB, IIIA and IIIB are known. We focus our attention on finding the metric in the 
region IVA {0 < u < L and < v < L), where the gravitational spin-spin interaction begins. 
If the value of j is small, we can expand the metric in terms of j. The background spacetime 
is the Minkowski spacetime, because the regions I, HA, and IIIA are fiat for j = 0. The 
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first-order perturbation is easily solved. Because the metric in the regions IIA and IIIA is 

{—dudv + dr"^ + r'^dcf^ + 2e{ulr)drd(f)^ region IIA, 
(59) 
—dudv + dr'^ + r'^dcf^ + 2e{av/r)drd(f), region IIIA, 

where 

e = 2j/L, (60) 
the metric in the region IVA is found to be 

ds^ = -dudv + dr^ + r'^dcj? + 2e drdcj), (61) 

r 

using the linearity of the first-order perturbation. Strictly speaking, we have to specify the 
properties of matter interaction between the sources of two incoming gyratons in order to 
determine the metric of the whole region IVA. The domain where the matter interaction 
is important is within the lightcone of the source collision, i.e. uv > r^. The spacetime 
structure of the other domain uv < in the region IVA is not affected by the matter 
interaction and therefore the metric (EH) can be applied for this domain. In the following, 
we restrict our attention to the domain uv < and do not consider the effect of matter 
interaction. 

In order to study the condition of AH formation, the above first-order solution is not 
sufficient because the nonexistence of the AH is due to the effect of nonlinear terms in j. 
Thus we should study (at least) the second-order perturbation, with which we will proceeded 
in this section. 

This study has the following meanings. First, it will clarify to what extent the condition 
of AH formation depends on the choice of the slice. In the previous section, we found that 
the conditions are very different in the (2a) and {2b) cases. Although we expected that this 
is due to the artificial effect of the slice choice, the study in this section will explicitly show 
whether such an expectation is correct or not. Next, by comparing the two cases a = ±1, 
we can study the properties of the gravitational field generated by the spin-spin interaction 
in the gyraton collision. As a result, we will find the dependence of the AH formation on the 
relative helicities of the incoming gyratons. For the old slice u > = v and v > = u, we 
found no difference between a = ±1 cases because the function G{u,r) on the chosen slice 
depends only on j^. However, the second-order structure of the region IVA will depend also 
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on a and it will lead to different conditions for the AH formation on the new slice, which 
consists of the future boundaries of regions IVA, IIB and IIIB as illustrated in Fig. [TH 
The gravitational spin-spin interactions were studied for a spinnin g te st particle around a 



rotating body 23|, for a massless particle passing by a rotating body 2J, |25|], and for binary 
systems of weakly gravitating bodies 3, 27|. In the case of binary systems, the contribution 
of spin-spin interaction to the relative acceleration, ass, between two bodies was calculated 

as 



oss 



n{Si ■ S2) + Si{n- S2) + S2{n ■ Si) - 5n{n ■ Si){n ■ S2) 



(62) 



where is a relative location f = fi — ^2, n := r/r, /i is the reduced mass. Si and S2 
are spins of two bodies. For Si = Sin and 5*2 = 5*2^?, the spin-spin acceleration becomes 
3ss = {Q/ ^r'^)SiS2'n. Therefore, for a binary with both spins aligned with the relative 
location vector r (i.e., both Si and S2 are positive), the spin-spin interaction is repulsive. 
Another example where the spin-spin interaction plays an important role is the Hawking 
emission of massless particles with spin (e.g., photons and neutrinos) by a Kerr black hole. 
In this process, the flux of particles with a given helicity created by the rotating black hole 
has anisotropic distribution. The black hole radiates more particles in the direction where 
the spin is aligned to the angular momentum of the black hole than in the opposite direction 



28, kiOll and see also 



3l| for higher- dimensional cases). This indicates the existence of 



a spin-spin interaction between the black hole and an emitted particle, which is repulsive 
when two spins have the same direction. If we assume that a spin-spin interaction similar to 
the above examples is present for a system of two relativistic spin particles, the black hole 
formation in the head-on collision of two gyratons with the same spin direction (cr = -|-1) 
is expected to be more difficult than that with the opposite spin directions (a = —1). The 
calculations in this section will confirm this. 

In Sec. V A, we derive the second-order Einstein equations and solve them. Then the AH 
equation and the boundary condition on the new slice is studied in Sec. V B. We present 
the numerical results in Sec. V C. In Sec. V D, we discuss the properties of the gravitational 
field in region IVA in more detail using the null geodesies. This helps us to interpret the 
results of AH formation. 
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A. Second-order equations 



We adopt e = 2{j/L) as a small expansion parameter and assume the following metric 
ansatz in {u,v,r,(j)) coordinates: 

11 ~\~ (Jl) 

ds^ = -(1 + e^c)dudv + (1 + e^a)dr'^ + r'^{l + e%)d(j)'^ + 2e— drd(j). (63) 

r 

Here a, b and c are functions of u, v and r. Expanding the Einstein equations up to the 
second order in e, we obtain^: 

/ V r / r 

Cnr + + ^ = ^ - + 0^^^), (67) 

c,vr + b^vr + — = — - ^(m + av), (68) 



.4a 

2Cr-T - ^a-uLV + b^rr + -&,r = — + —{u + (Tt;)''^ -, (69) 



2— - 46,„, + + -6,, = ^ + i.(u + avf - ^. (70) 



These relations follow from uu, vv, uv, ur, vr, rr, 00 components of the equation R^^ = 0, 
respectively. The other components vanish automatically. 

The initial conditions for this system are found by expanding the exact metric in regions 
IIA and III A in terms of j: 

a = b = c = 0, for = 0, (71) 

2 2 

a = b = c = 0, for u = 0. (72) 

The solutions satisfying these initial conditions are found as 

a = ^(n + avf + ^ (s " j^) + log(l - x)] , (73) 



^ More strictly, we should put grcj, = e[u9{u) + av9{v)]/r in Eq. ([63|) . although we do not show this because 
the equations become tedious. We note that this step function formula leads to the junction conditions 
atM = 0<w<L and v — < u < L through Ruu — Rw — 0, which the solution (|75)) - ([75)) satisfies. 
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b 



1 

4^ 



2r2 [ V 



1 — X 



1 



) 



+ log(l - x) 



(74) 



X 



(75) 



c = 



2^2 1 — 2;' 



with 



X = uv/r'^. 



(76) 



We discuss now the properties of the second-order solution in region IVA. First, a hne 
u,r,(f) = const is a null geodesic, although when c 7^ the coordinate v is no longer an affine 
parameter along the geodesic. Similarly a line v,r,(f) = const, is a null geodesic, although 
the coordinate u is not an affine parameter along it. Thus, the coordinates {u,v,r,(j)) 
simultaneously label the two null geodesic congruences. 

Next, all second-order quantities a, b, and c diverge at x = 1, i.e., uv = r^. Therefore it 
is interesting to ask whether x = 1 is a physical singularity or a coordinate singularity. For 
this purpose, we calculated the leading term in the expansion of the Kretchman invariant 
K := Rahcd.R"'^'^ near this point: 



Evidently it is divergent at x = 1. Because we are studying perturbation, the formula 
( 1771) cannot be trusted in the neighborhood of x = 1, and we cannot definitely claim that 
there is a physical singularity at x = 1. Still, Eq. (1771) indicates that there always exists 
the region where the perturbation breaks down around x = 1 for any small j. Hence, it is 
natural to expect that the exact solution, if it is found, also will have a real singularity of 
which location is shifted by O(j^) from x = 1. If this is the case, a physical singularity is 
produced atu = f = r = Oby the collision of gyratons and expands (almost) at the speed 
of light because uv = represents a light cone in the background spacetime. We note that 
this singularity formation is a consequence of the infinitely narrow shape of the source, i.e. 
Eqs. ([3]) and (jlj). In a realistic case where each source has a finite radius f = fg, the metric 
is regular at the source and therefore the singularity is not produced at uv = r^. Then 
the spacetime structure in the region uv > is determined by matter interaction between 
the two sources. Although the dependence on the properties of matter interaction is an 
interesting issue, it is not tractable by the perturbation. 

Finally, although the metric is continuous everywhere, ffist derivatives of a, b, and c are 
discontinuous atM = 0<f<L and v = < u < L. As a result, some components of 
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FIG. 15: The new slice on which we study the AH formation and the schematic shape of AH in 
the new shce. 



Riemann curvature, Rumr, Rvrvr, 



Ru 



RuAuA and RvSvA, have the delta function singularity there: 



¥4 



cr 



v'^{?,-x) 
r^(l — x) 
t;2(3-x) 
72(1 



X] 



9{u) - 2vS{u) 
e{u) + 2v6{u) 



e{v), 

e{v), 



(78) 
(79) 



and Rvrvr and Rv(f,v^ are obtained by changing u and v in Eqs. f l78|) and fl79l) . respectively 
(but note that Ricci tensor is zero in the sense of distribution). Hence, at the encounter of 
the two spin flows, a new shock field is produced and it grows linearly in u or v. The above 
four components of Riemann curvature are proportional to a and thus the feature of the 
shock gravitational field a.tu = 0<v<L and v = < u < L depends on the sign of a. 



B. AH equation on the new slice 

Because the metric in region IVA has properties that are somewhat different from other 
regions, we should derive the AH equation on the new slice. But the basic idea is the same 
as that in Sec. II. 

The second-order metric can be written like 

ds'^ = -Cdudv + Adr'^ + B{d4> + Ddrf, (80) 
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where A, B, C and D are functions of u, v and r. In region IVA, 



A = l + e 
and in region I IB, 



a — 



, B = r\l + e%), C=l + e^c, D = e 



[u + cri;) 



(81) 
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1, 



(82) 



where ul := u — L and we have kept terms up to second order in j. Based on this metric, we 
solve the AH equation on the new slice as shown in Fig. [151 The new slice consists of four 
parts: (1) = 0, L < u; {2) u = L, < v < L; (3) v = L, < u < L; (A) u = 0, L < v. 
Correspondingly, the AH consists of u = h^^\r) on the slice (1); v = h^'^\r) on the slice (2); 
u = h^^\r) on the slice (3); f = h^^\r) on the slice (4). 

We consider the AH equation for h{r) := h^'^\r). The tangent vector of the null geodesic 
congruence from the surface is given by 
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CD 
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and the condition of zero expansion becomes 
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or equivalently 
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(83) 



(84) 



0. (85) 



The equation for h^^\r) is obtained by just changing u and v in Eq. (!85|) . 

Now we explain the boundary conditions. At the intersection of the AH and the coordi- 
nate singularity m = L + r^, we impose 



(86) 



with the same reason as that in Sec. II. At the intersection of slices (1) and (2), we impose 
the condition that two null vectors of both sides of the surface be parallel, which is equivalent 
to 

h'^l^hf=4A, (87) 
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FIG. 16: The critical lines for AH formation on the (-L,j)-plane for the old slice. The results 
obtained by the exact formula and by the second-order approximation (denoted by "Exact" and 
"Approx.", respectively) are compared. The two results agree well and the error is 0{j^). 

where we used C = 1 on v = 0. Similarly, at the intersection of slices (2) and (3), we impose 

h^h^^^ = 4A/C. (88) 

In the cases cr = ±1, the functions A,B and C are symmetric with respect to u and v. 
Because D does not appear in the AH equation and the boundary conditions, the AH shape 
is symmetric with respect to the plane u = v. Hence, we only have to study h^^\r) and 
h^'^\r), and the boundary condition (|88|) is reduced to 

h[^^ = 2,/A/C. (89) 

We also note that because the functions A, B, and C do not depend on the sign of j, the 
condition for the AH formation is written in terms of |j| and L for each a. For this reason, 
j is assumed to be positive without loss of generality in the following. 

The numerical procedure is as follows. First we choose some value of tq and start solving 
/i(i)(r) with the boundary condition h^^\ro) = L + Tq and ( 186|) . When h^^\r) becomes L 
at r = ri2, we solve /i'-^-* using the boundary conditions h'^'^\ri2) = and h^? = 4A/h^r^. 
When h^'^^ becomes L at r = r23, we check whether the boundary condition (!89l) is satisfied. 
Iterating these steps for various values of ro, we can judge the existence of the AH and find 
its location. 

C. Numerical results 

In order to test the reliability of the second-order approximation, we studied the condition 
of AH existence on the old slice u > = v and v > = u using the exact formula and 
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FIG. 17: The top view of the AH (sohd lines) and the inner boundary of the trapped region (dashed 
Unes) in the case a = +1 for L = 0.3 and j = 0.05, 0.1, and 0.1129. 
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FIG. 18: The top view of the AH (solid hnes) and the inner boundary of trapped region (dashed 
lines) in the case a = —1 for L = 0.3 and j = 0.1,0.15, and 0.1581. For j = 0.1, there is no inner 
boundary of trapped region. 

the second-order formula for G{u,r) and compared the two results. Figure fT6] shows the 
regions of AH formation on the (L, j)-plane obtained by the two methods. The two results 
agree well and the difference is O(j^). This demonstrates the reliability of the second-order 
approximation for the old slice. Later, we will discuss also the reliability of the approximation 
for the AH study on the new slice. 

Now we show the results of the new slice. We first show the case a = +1 where the 
spins of two gyratons have the same direction (i.e., the helicity of one gyraton is positive 
and that of the other is negative). Figure [T71 shows top views of the AH shape for L = 0.3 
and j = 0.05, 0.1 and 0.1129. We could not find the solution for j > 0.1130. Similarly to the 
case of the old slice, there is some critical value of the spin j^'^:^l{L) for the AH formation. 
There are two solutions for each j < j^'^^llL), which correspond to the AH and the inner 
boundary of the trapped region, and they coincide at j = j^^^ (L). 

Next we show the case a = —1, where the spins of two gyratons are oppositely directed 
(i.e., the helicities of gyratons are both positive or negative). Figure fT8l shows top- views of 
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FIG. 19: The critical lines for the AH formation for the new slice for a = +1 (lower solid line) 
and for cr = —1 (upper solid line). The critical line for the old slice is also shown (dashed line). In 
both cases a = ±1, the AH formation is allowed in a larger region on the (L,_7)-plane compared to 
the old slice. The allowed region of cr = —1 is larger than that of cr = +1. In the case cr = — 1, the 
perturbative quantity becomes large at L > 0.9 and the shape of the critical line cannot be trusted 
there (dotted line). 

the AH shape for L = 0.3 and j = 0.1, 0.15, 0.1581. For j < 0.1, only one solution is found. 
Hence there is an AH but no inner boundary of the trapped region. For 0.15 ^ J < 0.1581, 
two solutions are found for each j. Thus the inner boundary of the trapped region appears 
for these values of j. For j > 0.1582, there was no solution. The critical value j^^l{L) of AH 
formation for cr = — 1 is larger than icrit (-^)- Thus in the case cr = — 1, the AH is allowed to 
form in a larger parameter regime compared to the case cr +1. 

Wc studied the critical value j^^j^(L) as functions of L for the cases c = ±1. Before 
showing the obtained results, we comment on the reliability of the second-order approxi- 
mation. In order to evaluate the error, we checked the maximum values of e^a, e^b and e^c 
on the AH at the critical line. In the case cr = +1, max[e^a] < 0.32, max[e^6] < 0.20, and 
max[e^c] < 0.14 are satisfied for arbitrary L. Therefore the expected error is about 20%. In 
the case cr = —1, they are found to be max[e^a] ~ 0.38, max[e^6] ~ 0.13 and max[e^c] ~ 0.12 
for L < 0.3, and thus the expected error is about 30% for L < 0.3. However, for L = 1.0, 
their values grow large: max[e^a] ~ 0.96, max[e^6] ~ 0.97 and max[e^c] ~ 0.95, and the 
approximation obviously breaks down at L > 0.9. Thus unfortunately we cannot trust the 
shape of the critical hne for L > 0.9 in the case cr = — 1. To summarize, we can trust the 
shape of the AH critical line of cr = -|-1 for arbitrary L and that of cr = —1 for L < 0.3 with 
the error discussed above. 
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Figure [T9l shows the parameter regions in the (L, j)-plane that allows the AH formation 
on the new slice in the cases cr = ±1, together with that on the old slice. In both cases, 
icrit(-^) go6s to zero in the limit L ^ 0. The critical line of a = +1 crosses the L-axis at 
L = 1, which is much larger compared to L = 1/e in the case of the old slice. Although the 
error in j^~J^(L) grows large for L > 0.9, the critical line of a = —1 does not seem to cross 
the L-axis for < L < 1. Hence, the allowed regions of the new slice is much larger than 
that of the old slice for both a = ±1. At the end of the previous section, we stated our 
expectation that the large difference between the allowed regions of the (2 a) and (26) cases 
is due to the artificial effect of the slice choice. It is now confirmed, since the allowed region 
of the (26) case has become much larger by just changing the slice. Comparing the two cases 
a = ±1, Jcrit (i^) is greater than jctili-^)- Therefore, the AH formation in the case a = — 1 is 
allowed in a larger parameter region compared to the case cr = +1, and the condition of the 
AH formation depends on the relative helicities of incoming gyratons. To briefly summarize, 
on the new slice, the condition of the AH formation is roughly expressed as L ~ 0.5 and 
7 < 0.1 in the case a = +1 and L ~ 0.5 and j < 0.15 in the case a = —1. 

The reason why the allowed region in the case cr = +1 is limited to < L < 1 is 
understood as follows. In the case j = 0, the AH solution is given by 

h'-^\r) = h'-^\r) = L + 2r^ log {r/r 12) , (ri2 < r < v^ria), (90) 

/i(2)(r) = h^^\r) = 2 log (r/ria) , {r^ < r < 1), (91) 

with = 1/ \/L. Although the AH is expected to converge to this solution in the limit 
j — ^ 0, we should take care of the presence of the singularity uv = r"^, where the perturbative 
quantities diverge. For L > 1, the singularity crosses the surface fl9T|) . invalidating it to be 
an AH. Hence j = and j = 0^ are different for L > 1, and no AH exists for small j. On 
the other hand, for < L < 1, the surface fl^U]) and fPT]) is the AH in the limit j ^ 0, 
because the singularity does not cross the surface. Hence, it is natural that the region of 
the AH formation is restricted to < L < 1. Although we could not specify the allowed 
region for a = — 1 around L ~ 1, the above discussion would hold also for this case. Hence, 
if the exact solution for region IVA is found, the allowed region for a = —1 will turn out to 
be restricted to < L < 1 

^ This discussion holds only for a collision of gyratons with singular sources, Eqs. ^ and In the collision 



35 



D. Gravitational field in the region IVA 



We discuss the properties of gravitational field in region IVA in more detail, because it 
helps us to understand the reason for the different allowed regions in the cases cr = ±1. For 
this purpose, we study the "gravitational force" acting on the null geodesies m, r, = const, 
and v,r,(f) = const. 

Let us consider a null geodesic congruence u = Uq, r = tq, < (f) < 1. The section of the 
congruence and v = const, is a loop and the quantity 

rl2;:'\v)=ro{l+e\uo,v,ro)/2) (92) 

gives a radius of the loop (i.e., the proper circumference divided by 2tt). We define the 
"gravitational force" F^'^°'^'^\v) toward the symmetry axis by 



n2„("o,r-o) 



^K,ro)(^) := . (93) 



The force is attractive if F^'^°'^°\v) < and repulsive if F^'^°'^°\v) > 0. Similarly we consider 
another congruence v = vq, r = r^, < (f) < 1 and introduce its loop radius 'rl^op°\'^)- Then 
another kind of force is defined by 

F(^0'^o)(m) := . (94) 



The two forces are calculated as 

2auS{v) + 1 + 0- 



2av5{u) + 1 + a- 



(95) 



(96) 



r^(l — xy 

The delta function of the first term in the square brackets of each formula comes from the 
new shock field at-u = 0<t'<L and v = < u < L [see Eqs. ( 1781) and (1791)]. 

In the case a = +1, both F^'^''^\v) and F^'"''^\u) are positive outside of the singularity 
X = 1. Hence, the gravitational field is repulsive in the whole region IVA. On the other 
hand, in the case a = — 1, the coefficients of the delta functions in Eqs. (195|) and (196|) are 
negative, indicating that the new shock fields are attractive. The third term in the square 



of realistic beam pulses, the singularity is not produced at uv = and the regions of AH formation might 
become different from Fig. ITHl 
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FIG. 20: The behavior of the function b{uQ, v, tq) in the cases a = +1 (left) and a = —1 (right) for 
ro = 1 and uq = 0,0.5, 1. The value of b{uo,v,rQ) is directly related to the radius '^ioop'''^^('^) °^ 
light ray congruence u = uq, r = tq, < (p < 2tt. The light rays quickly bend at v = due to the 
delta functions in Riemann curvature (|78|) and (|79|) in both cases but the bending directions are 
opposite. b{uo,v,rQ) continues to increase in the case a = +1, while its behavior strongly depends 
on uq in the case a = —1. 

brackets is also negative. If x is close enough to 1, the third term exceeds the second term 
and the force becomes negative. Hence, around the singularity x = 1, there is always the 
attractive region of the gravitational force. If x is close to 0, the third term is smaller than 
1 and the gravitational field is repulsive in such a region. Therefore, both attractive and 
repulsive regions exist for a = —1. 

Let us look at the behavior of the loop radius r[^°p°^ (v) . Ignoring a factor, the change in 
r[^op°\'^) is presented by b{uo,v,ro). Figure [20] shows the behavior of b{uo,v,ro) for tq = 1 
and Mo = 0, 0.5, 1 for the two cases a = ±1. Because of the delta function in the force (p5ll . 
b{uo, V, ro) is not smooth at t> = for no > in both cases. In the case a = +1, b{uo, v, ro) 
suddenly increases at f = and blows up, since the force is repulsive everywhere. In the 
case cr = —1, b{uo,v,ro) suddenly decreases at f = due to the attractive force. For t> > 0, 
the behavior of b{uQ, v, ro) strongly depends on the value of uq- If uq is large, the light ray 
feels repulsive force at the beginning but later feels attractive force. 

Now we discuss the reason for j[J^{L) > j^^:^l{L), i.e., the difference between the allowed 
regions for the AH formation of o" = ±1. In the case a = +1, the gravitational field in 
the region IVA is repulsive everywhere. If j is increased, the repulsive force exceeds the 
attractive force generated by the energy, causing the disappearance of the trapped region. 
On the other hand, in the case a = —1 there are both attractive and repulsive regions. 
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FIG. 21: The sign of two forces and on slice (2), i.e. < i; < L = it, for 

a = —1. The slice is divided into three regions: region ( ) where F^'^''^^ (v) < F^'"''^\u) < 0, 

region (-+) where ^("'^^(v) < < F'^'"'''\u), and region (++) where < < F(^'^)('u). 

The unit of r (horizontal line) and v (vertical line) is the spin duration L and the gray line indicates 
the singularity x = 1. 

Figure EH shows the sign of two forces F^'^'^\v) and F^'"'^\u) on slice (2), i.e. < v < L = u. 

The slice is divided into three regions: region ( ) where F^'^'^\v) < F^'"'^\u) < 0, region 

(-+) where F("''')(w) < < F(^''')(m), and region (++) where < F("''')(t;) < ^(''■^^(u). 
This figure shows that the gravitational force is attractive around the singularity x = 1 
and repulsive for r ^ L. For L <^ 1, the attractive region is a tiny portion just around 
the singularity and the force is repulsive almost everywhere on the surface. The repulsive 
force becomes dominant as j is increased, resulting in disappearance of the AH. Although 
the attractive region becomes large for L ^ 1, the attractive force does not help the AH 
formation effectively since there is the constraint Jcrit(-^ = 1) = coming from the size of 
the singularity a; = 1 as discussed in the previous subsection. Therefore, also in the case 
a = —1, the spin j makes the AH formation more difficult. However, in the case a = —1, 
the repulsive force is obviously smaller than that of the case a = +1 for a fixed j value. 
Hence a larger value of j is needed for the disappearance of the AH. This explains our result 

VI. SUMMARY AND DISCUSSION 

In this paper, we studied the AH formation in the head-on collision of gyratons. We 
introduced four gyraton models in Sec. II: a spinless p-gyraton, an AS-gyraton, and spinning 
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TABLE I: Summary of the obtained results. For each case, the condition of AH formation was 
found in terms of L [the energy duration in the (0) case and the spin duration in other cases] and 
the spin value j = J/2prh{2p) (assumed to be positive). The unit of the length is rh{2p) = AGp. 



collision type 


slice ((j) 


gyraton 1 


gyraton 2 


condition 


of AH formation 


(0) 




p 


p 




L<1A 


(la) 




a 


AS 


L ~ 1 


j < 0.25 


(16) 




b 


AS 


L ~ 1 


j < 0.25 


(2a) 




a 


a 


L ~ 1 


j < 0.2 


(26) 


Old 


b 


b 


L ~ 0.15 


3 < 0.1 


(26) 


New (+1) 


b 


b 


L ~ 0.5 


3 < 0.1 


(26) 


New (-1) 


b 


b 


L ~ 0.5 


3 < 0.15 



a- and b- gyratons. The energy and spin profiles of each gyraton are given in Eqs. fl71)- f|T0l) . 
For a spinless p-gyraton and an AS-gyraton, the energy profile is a step function with width 
L and a delta function, respectively. For a- and b- gyratons, the energy profile is a delta 
function and the spin profile is a step function with width L. The difference between a- and 
b- gyratons is the relative locations of the energy and spin profiles. We introduced the null 
geodesic coordinates for each gyraton, and discussed the property of its gravitational field. 
Especially, a spinning gyraton has a repulsive gravitational field around its center. 

Then the problem of the head-on collisions of two gyratons was set up and the AH was 
studied on the slice u = > v and t> = > m in Sees. II-IV. The studied collision cases 
and obtained results are summarized in Table [B In all cases two gyratons are assumed to 
have the same energy. Case (0) is the collision of two identical spinless p-gyratons. In this 
case the energy duration L should be smaller than some critical value for the AH formation. 
In cases (la) and {lb), we studied the collision of spinning a- and b- gyratons with an 
AS-gyraton, respectively. We obtained the conditions for the AH formation in terms of the 
spin duration L and the spin j. They are shown in Fig. [10] and roughly summarized as 
in Table [B (Here j > is assumed since the AH formation does not depend on the spin 
direction.) In both cases, there was a critical value jcrit(-^) for the AH formation for a given 
L. We found no significant difference between the two cases. In cases (2a) and {2b), we 
studied the collision of two spinning a- and b- gyratons, respectively. Two gyratons were 
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assumed to have the same spin duration L and absolute value of the spin j. We obtained 
the conditions for AH formation in terms of L and j. They are shown in Fig. [13] and roughly 
summarized as in Table [H (Here j > is assumed and the relative direction of two spins a 
is not specified, since the AH formation does not depend on the directions of two spins on 
the studied slice.) We found that the allowed region on the (L,j)-plane in the {2b) case is 
much smaller than that in the (2 a) case. 

In Sec. V, we focused our attention on the gravitational spin-spin interaction after colli- 
sion in the (26) case. We solved a part of the future to the slice u = Q < v and f = < m 
(old slice) in the collision of gyratons with spins j and aj (a = ±1), using a method of 
perturbation where j is a small expansion parameter. The solved region is the past to the 
collision of the energy flows, but the two spin flows interact with each other in that region 
(see Fig. [H] for details). Therefore we could study the spin-spin interaction. Then we again 
studied the AH formation on the future edge of the solved region (new slice) and compared 
the obtained results to those of the old slice. It was found that the allowed region becomes 
larger by just changing the slice (Fig. [T9ll . Hence, the difference between the results of old 
slice in cases (2 a) and (26) was due to the artificial effect of choosing a slice on which we 
study the AH. 

Furthermore, we found the dependence on the relative helicities of incoming gyratons. In 
the case a = +1 where two spins have the same direction (i.e., helicities have opposite signs), 
the gravitational field is repulsive everywhere due to the spin-spin interaction. On the other 
hand, in the case a = —1 where two spins have the opposite directions (i.e., helicities have 
the same sign), the spin-spin interaction decreases the repulsive force and even changes it 
into the attractive force in a part of the studied region. Correspondingly the allowed region 
of the AH formation for a = — 1 is larger than that for o" = +1 (Fig. [T9l) . 

In the light of the above studies, we claim the following. For the AH formation in the 
head-on collision of gyratons, (i) the energy duration should be smaller than some critical 
value (close to the system gravitational radius); (ii) the spin duration should be at least 
of order of the system gravitational radius (it should not be too small or too large); (iii) 
the spin should be smaller than some critical value that is a function of the spin duration. 
Further, (iv) the AH formation in the collision of two gyratons with the oppositely directed 
spins is easier than that with the same direction of spins. 

Now we discuss the possible applications of the obtained results for mini black hole 
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production at the LHC in the TeV gravity scenarios where Mp = TeV. Let us consider 
the colhsion of two spinning particles, and use our result of the (2a) case, i.e. the collision 
of two identical a-gyratons, for the condition for the black hole formation as an example. 
Restoring the length unit, it is written as L ^ rh{2p) and J < 0.4 x p rh{2p). We use the 
Lorentz contracted proton size L ~ 1.5 x lO^^fm for the spin duration and put J = h/2 as 
possible candidates for these values. Substituting p = (few)Mp and rh{2p) = {fevf)h/Mp, we 
find L ~ rh{2p) and 

0.4 X p rh{2p) ~ (few);i > h/2 = J. (97) 

Hence, the above two conditions are satisfied and the black hole is expected to form in the 
head-on collision under our assumption. Thus, the effect of spins of incoming particles might 
not be significant for the black hole production rate. Still, the spin might change the cross 
section of the black hole production by a factor and studying this effect would be interesting. 



We also revisit the study by Giddings and Rychkov [ll|, because our study is related to 
the assumption they made. In that paper the collision of quantum wave packets with width L 
was considered. Their result is that if h"^/ {rhM^) <^ L <^ rh-, the higher- curvature correction 
is small and the predictions by general relativity are reliable. The latter inequality L 
was imposed by the expectation that the gravitational field of such a wave packet would be 
sufficiently close to that of the AS particle and thus the AH would form in a collision of such 
two wave packets. Our result of the (0) case, i.e. the collision of two identical p-gyratons, 
explicitly demonstrates the accuracy of this expectation. Moreover, because we found the 
AH also for L < lAvh, the condition can be relaxed to fi^ / {r^M'^) <^ L < rh- Note that 
this criterion holds also for wave packets of spinning particles, if their energy is sufficiently 
large, p ^ Mp. Our results of the (2a) and {2b) cases, i.e. the collisions of two identical 
spinning a- and b- gyratons, show that ^ is necessary for the AH formation for small 
L. Restoring the length unit and adopting J = h/2, it is rewritten as ^^/(IGr/^p^) < L. 
However this does not provide a new condition since the original condition h"^/ {r^Mp) ^ L 
implies h? / [IQrhP^) < L for p ^ Mp. Therefore, our results do not contradict the claims in 



The important remaining problems are as follows. The first one is to explore the case 
a = —1 further. This is because the condition of the black hole formation is expected to be 
different from that of the AH formation. In the case cr = +1, however, the critical value of j 
for the black hole formation will remain finite, because both the gravitational field generated 
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by the spin source and the spin-spin interaction are repulsive. On the other hand, in the case 
a = —1, the repulsive gravitational field of each incoming gyraton is weakened and becomes 
even attractive in some part of the spacetime by the spin-spin interaction as shown in Sec. 
IV. Hence, there is the possibility that later the gravitational field turns to be attractive 
everywhere and the critical value of j blows up. 

The next problem is the collision of gyratons with a nonzero impact parameter. In 
these grazing collisions, new effects of the spin-orbit interaction will appear. Moreover, the 
properties of spin-spin interaction might change. Let us recall Eq. (!62|) . the acceleration 
ass due to the spin-spin interaction between weakly gravitating bodies. In the grazing 
collisions, the spins are orthogonal to the relative location vector and ass is calculated as 
ass = ~(3//i?"^)(>5i ■ S2)n. Therefore in the aligned (resp. anti-aligned) case, the spin- 
spin interaction becomes attractive (resp. repulsive), which is opposite to the head-on 
collision case. Therefore it is expected that the nonzero impact parameter would make the 
interactions more complicated but more interesting. 

It is also important to simulate the collision of gyratons with realistic sources. In this 
paper, we assumed that each incoming gyraton has a singular source, Eqs. ([3]) and (jll), 
and studied only the spacetime regime where the matter interaction is not important (i.e., 
uv < in Sec. V). In a realistic situation, however, the source of an incoming gyraton is a 
beam pulse with a finite radius fg. Then, the matter interaction determines the spacetime 
structure within the lightcone of the source collision, and the condition for the black hole 
formation will depend on the properties of matter interaction. In order to study this effect, 
we should solve the Einstein equations together with the field equations for the sources. 

Finally, the generalization for the higher- dimensional case is necessary to obtain the 
results that can be directly applied for the black hole production at accelerators in the TeV 
gravity scenarios. 
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